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Abstract 

The study of the response of the two modules of the large-angle electromagnetic 
shower calorimeter (LAC) of the CLAS detector to charged and neutral particles 
is reported. The results agree very well with the Monte Carlo simulation. The pro- 
cedures adopted for the energy and timing calibration are also discussed, proving 
that the module geometry allows for simple self-calibrating energy and timing al- 
gorithms. 
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1 Introduction 



The accelerator CEBAF at Jefferson Lab provides up to 6 GeV electrons to 
three experimental halls. Hall B houses the Large Acceptance Spectrometer 
CLAS, which is designed to study multi-particle final-state reactions induced 
by photons and electrons at luminosities up to 10 34 cm~ 2 s _1 . The charged 
particle momentum is analyzed using a 1 T toroidal magnetic field with az- 
imuthal symmetry generated by six superconducting coils arranged around 
the beamline. Each sector between two coils acts as a single spectrometer 
instrumented with three layers of drift chambers (DC) for track reconstruc- 
tion, one layer of scintillator counters for time-of-flight measurement (TOF), a 
gas filled Cherenkov counter (CC) and an electromagnetic shower calorimeter 
(EC) . Drift chambers and time-of-flight define the sensitive region of each sec- 
tor in the range of the polar angle from 10° to 150°, the forward calorimeter 
modules cover from 10° to 45° in each sector and the two modules of the large- 
angle calorimeter (LAC) extend the coverage up to 75° in two CLAS sectors. 
In 1998 the two LAC modules as a part of the CLAS detector have been com- 
missioned and since then they have been successfully operating to run the Hall 
B experimental program with electron and photon beams. The main features 
of the LAC module have been already described in previous papers where we 
reported the properties of the main components as the plastic scintillators 
[2,4], the light guides [3] and the photomultipliers [5]. The detector structure 
and the response to minimum ionizing particle has been described in Ref. [6]. 
In this paper we report on the response of the LAC to charged and neutral 
particles. The results are based on data taken with both an electron and a 
photon beam impinging on a cryogenic target filled with liquid Hydrogen or 
liquid Deuterium. We first discuss the energy and timing response obtained 
from the detection of pions and protons. After that we show the performances 
for neutral particles: photons and neutrons. 



2 The LAC Energy Response 

To discuss the energy and timing responses, we first give a brief description 
of the detector. The conceptual drawing of the internal structure is shown 
in Fig. 1: the LAC module has a rectangular shape with a sensitive area of 
217x400 cm 2 and the internal structure is made of 33 layers, each composed 
of a 0.20 cm thick lead foil and 1.5 cm thick NE110A plastic scintillator bars. 
Each scintillator layer is protected from the contact with the lead by thin 
Teflon foils. The width of the scintillators is roughly 10 cm and is slightly 
increasing from the inner layers toward the outer layers to provide a focusing 
geometry. Scintillators in consecutive layers are rotated by 90° to form a 40x24 
matrix of cells with area approximately 10x10 cm 2 . The module is vertically 
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Fig. 1. The conceptual drawing of the internal structure of the LAC module. 
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Fig. 2. The cell selection procedure: minimum ionizing particles release energy in a 
region that is generally inside the dark grey square defined by the crossing of two 
orthogonal stacks (shown in light grey). The emitted light propagates along each 
stack and is collected at both end of the two stacks generating four ADC signals. 



divided into an inner (17 layers) and an outer (16 layers) parts with separate 
light readouts. Scintillators lying one on top of the other with the same orien- 
tation form a stack. The focusing geometry is such that any straight trajectory 
that originates from the target would hit two crossing stacks both in the inner 
part as well in the outer part. For each stack the light is collected at both 
ends separately using light guides coupled to EMI 9954A photomultipliers. 
For each module there are 128 stacks and 256 photomultipliers. 
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2.1 The Energy Calibration 



The energy calibration is performed continuosly checking off-line each pho- 
tomultiplier. The ADC pedestals are measured every day and the gain in- 
homogenities are corrected off-line setting the peak of the minimum ionizing 
particle (MIPs) distribution to the same value for all photomultipliers accord- 
ing to the following procedure already used during the tests with cosmic rays 
[6]. Minimum ionizing particles have enough kinetic energy to penetrate in 
the detector and exit from the back side. The events corresponding to those 
particles (mainly high-momentum pions) that release their energy in a single 
cell are selected, in the off-line analysis, by imposing the presence of four ADC 
signals in the inner part and the four corresponding ADC signals in the outer 
part as shown in Fig. 2. As a preliminary step, the correction for the depen- 
dence from the cell position is applied considering that the light attenuation in 
scintillators relates the read-out charge Q(x) to the "true" value Q according 
to the expression: 

Q(x) = Q (e' x ' x + ae-( 2L -*)/ A ) , (1) 

where x is given by the coordinates of the center of the cell, L is the stack 
length and the two parameters A and a (representing, respectively, the light 
attenuation length and the reflectivity of the end surface of the scintillator 
bars) are the constants that have been determined for each stack by the cos- 
mic ray calibration described in Ref. [6] and stored in a database. After the 
position correction, the conversion factor charge-to-energy is evaluated, for 
each photomultiplier, in order to get the same MIPs peak position and to 
minimize the width of the overall deposited energy distribution. 

After the energy calibration and the equalization of the read-out of the 512 
photomultipliers, the calorimeter energy response has been studied for differ- 
ent detected particles and the results are discussed in the following sections. 

2.2 The Response to the Ionizing Particles 

Protons are a powerful tool to check the linear response of individual stacks in a 
wide energy interval. The protons detected by CLAS and used in this analysis 
originated from multipion emission off a liquid hydrogen target in electron 
scattering. They have been selected by (3. vs. momentum discrimination using 
informations from TOF and DC. Fig. 3 shows a typical distribution of the 
deposited energy in the LAC as a function of the proton momentum. Protons 
with momentum less than 0.6 GeV/c release all their kinetic energy in the TOF 
placed in front of the LAC scintillators as well as in the LAC front-housing 
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Fig. 3. The energy deposited in the LAC module by protons .vs. the particle mo- 
mentum P. 

and are not detected in LAC. Protons with momentum up to 1 GeV/c are 
completely stopped by the calorimeter while higher-momentum protons release 
a fraction of their kinetic energy that approaches at very high momenta the 
asymptotic value of 106 MeV. This value corresponds to the energy released 
by MIPs in the LAC. 

The momentum dependence of the deposited energy and the width of the de- 
posited energy distribution agree with the Monte Carlo simulations. In Fig. 4 
the momentum dependence of the width a E of the deposited energy distribu- 
tion of the protons stopped in the calorimeter is reported. The dependence of 
cte from the deposited energy E has been fitted with the expression 



with E in GeV . The parameter a = (2.5 ±0.1) * 10~ 3 GeV resulted very 
small while the value b = (18.4 ± 1.0) * 10~ 3 GeV 1 / 2 is larger than expected 
assuming the contribution only from the fluctuations in the light transmis- 
sion and collection. Actually different effects contribute to the widening the 
deposited energy distribution; they include non uniformity in the effective scin- 
tillators and lead thicknesses due to different proton incident angle as well as 
non-uniformity in scintillator intercalibration. The contribution of each single 
effect cannot be easily disentangled, therefore we compared the overall re- 



a E = a + bE 1/2 , 



(2) 
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Fig. 4. The deposited energy resolution vs. the particle momentum P for the protons 
completely stopped in the calorimeter. The continuous line represents the result of 
the Monte Carlo simulation. 

sponse of the LAC modules with a realistic Monte Carlo simulation where we 
assumed 5 photoelectron/ MeV collected in average by each photomultiplier 
as already measured in cosmic ray tests [6], and a photomultiplier intercali- 
bration non- uniformity equal to 5 %. The continuous line in Fig. 4 represents 
the result of the simulation and the agreement is fairly good. 



2.3 The Response to the Electromagnetic Showers 

To study the response to electromagnetic showers, the only possibility is the 
detection of the photons from the 7r° decay, because the flux of electrons 
scattered at large angles is too low and the 7r~ contamination too high. The 
photons are identified selecting those events where the LAC energy signal 
has no TOF and DC signals in coincidence. The size of the electromagnetic 
shower is normally larger than the stack width and several stacks are involved 
in the detection. To reconstruct the photon energy first the two coordinates 
of the center of the shower are determined by weighting the coordinates of 
contiguous non-zero ADCs with their read-out, then for each ADC the read- 
out is corrected according to Eq. 1. Finally the photon energy is obtained 
multiplying the sum of the corrected ADC read-outs by the sampling factor 
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/ = 0.32. This factor has been determine! from Monte Carlo simulations. 



Unfortunately the photon energy spectra shows a continuous yield and cannot 
provide a direct test of the energy calibration. Nevertheless photon energy res- 
olution could strongly affect the reconstructed mass of 7r° and its investigation 
can provide an indirect check of this procedure as discussed in Sec. 4. 



3 The LAC Time Response 

The time response of the detector can affects the detection and identification 
of neutral particles. Photons originated from mesons decay can be disentan- 
gled from neutrons only using their different time-of-flight. To provide 7 — n 
discrimination good timing uniformity and resolution are required. The de- 
tector geometry allows a self-calibrating procedure. For each photomultiplier 
the relation between the particle time-of-flight t and the corresponding TDC 
read-out t' already corrected for the time-walk according to Ref. [6] can be 
simply written as 



being / the length of the trajectory from the target to the detector, (3c the 
particle velocity, x the distance between the impact point in the detector and 
the photomultiplier, v sc the propagation velocity of the light in the scintillator, 
r the intrinsic time delay of each electronic chain and To the overall calibration 
constant equal for all photomultipliers. To provide the correct relationship 
between i! and f3, the 512 different r constants and T need to be set properly. 

We first proceeded to the equalisation of the time response by setting the r 
values using signal from MIPs. Only events hitting a single cell were selected 
imposing the presence of four ADC signals in the inner part and of the four 
corresponding ADC signals in the outer part of the module. The ADC sig- 
nals were originated from pairs of photomultipliers coupled respectively to the 
right and to the left hand side of the same scintillator stack. By using the 
corresponding TDC readouts t' Le f t and t' Right together with Eq. 3 we evaluated 
for each stack the half-sum 



t = l/p c = t'-r- x/v sc - T , 



(3) 



St — - (t mght + tLeft) 
_ ^Right "j~ txeft 



r Right + r Left 

2 



L/2v sc - T c 



(4) 
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being L the scintillator lenght, and the half- difference 
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Fig. 5. The At timing calibration scheme. For a couple of crossing stacks the coor- 
dinates and the calibrated value of At for the first cell are reported. The picture 
shows that with the assumed calibration, the At readout for the cell located at the 
crossing point is the same for the two stacks. 



At — - (t Right — t Le f t ) 
t' — ¥ 

_ L Right L Left 



(5) 



T Right ~ T Left 



With the adopted selection each event is represented by eight numbers: the 
values Sts and At s for the inner-short stacks, St^ and At L for the inner-long 
stacks and similar values for the outer stacks. St is simply the particle time- 
of-flight and consequently, for each single event, the four St have to assume 
the same value, or equivalently, to satisfy the equation 



St a — s^ — 0, 



(6) 



for any pair (a, b) of the St variables. 

A similar relationship can be obtained for the four At values as described in 
Fig. 5: for each stack the At value has been set such that for events in the 
first cell the correspondig timing would be 

At(0) = x/v sc , (7) 



being x the stack distance from the bottom-left angle of the module or equiv- 
alently the stack coordinate. In this way it is easy to verify that the events 
detected in any cell with coordinates (x, y) are associated to the same At value 
for both short and long read-out: 



x + y 

At s = At L = y -. 

Vsc 



(8) 
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Fig. 6. The distribution of the difference Sis— Siz, in the time-of-fiight reconstructed 
with inner-short and inner-long stacks of the LAC for charged pions and protons 
detected in both modules. 

With this condition, for each event the four At have to assume the same value, 
or equivalently 

At a - At b = 0, (9) 
for any pair (a, b) of the At variables. 

Therefore it can be easily proved that equalizing the timing response of the 
LAC is equivalent to finding the values of the 512 r that simultaneously set 
to zero the average value of the distributions (Etg — Et^) and (Atg — At^) 
separately for the inner and the outer part of the detector, minimizing at the 
same time their widths. In Fig. 6 a typical result is reported. 

To fix T we used the time-of-fligth of the photons originated from 7r° decay. We 
selected photons candidates in the reaction ep — > e'pX with missing mass equal 
to M w o and a neutral particle detected in coincidence by the LAC modules, 
then we evaluated from the time-of-flight the corresponding (5 = IjcLt values 
being / the length of the straight trajectory from the target center to the 
reconstructed impact point in the LAC. We set the T value imposing the 
mean values of the (5 distribution equal to 1. In Fig. 7 the distribution of (5 for 
detected photons shows a gaussian shape with 073 = 2.6%. Considering that 
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Fig. 7. The reconstructed (3 for photons detected in the LAC modules. 
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Fig. 8. The time-of-fhght reconstructed with TOF vs. LAC for charged pions and 
protons. 
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the average distance from the target is 4 m this value corresponds to a timing 
resolution o t = 250 ps well in agreement with that obtained in the preliminary 
tests with cosmic rays [6]. The timing response is very good in the whole (3 
range as shown in Fig. 8 where the (3 reconstructed with LAC is compared to 
the 13 from the TOF for protons {(3 < 0.85) and tt+ {(3 « 1). The data were 
taken using a 4.0 GeV electron beam incident on a liquid H 2 target. 



4 The Reconstruction of the 7r° — >• 27 Decay 

The reconstructed value of the 7r° mass provides a quantitative check of the 
LAC energy equalisation and calibration. Any energy miscalibration would in 
fact provide an incorrect mean value and a larger width for the reconstructed 
mass distribution. To verify that the energy calibration procedure discussed 
in Sec. 2 is correct, we selected the events from the reaction ep — > e'^X with 
one photon detected in the LAC modules and one in the EC modules. The 
neutral hits in the calorimeters were identified by the absence of a matched 
track in drift chambers. Photons detected by the LAC generally belongs to low 
momentum pions that decay emitting the two photons with a large relative 
angle 12 ; therefore candidate events corresponding to the decay it — > 27 were 
required to have a neutral hit in the LAC and a second neutral hit in the other 
LAC module or more favourably, in any EC module. 

The meson mass has been reconstructed from the photon energies Ei, E 2 and 
the angle 9 12 between the direction of the two photons according to 

M 7 2 7 = 2E 1 E 2 (1 - cos9 12 ) . (10) 

The reconstructed M v o, as shown in Fig. 9, peaks at the right energy and the 
width of the distribution a = 28% is well in agreement with the Monte Carlo 
simulations. The data were taken using a 4.0 GeV electron beam incident on 
a liquid H 2 target. 



5 The Neutron Detection 

The LAC detector has been designed to provide high neutron detection effi- 
ciency From Monte Carlo simulations a detection efficiency larger than 20% 
for neutron momenta higher than 0.7 GeV/c is expected. Neutrons as well as 
photons are identified by a neutral hit in the LAC and the particle velocity 
is reconstructed from the time-of-flight. Neutrons are separated from photons 
setting (3 < 0.95. 
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Fig. 9. Reconstructed invariant mass for photons from the reaction ep —> e'jyX 
with at least one photon detected in the LAC modules. 




0.3 0.4 

Fig. 10. The relative difference between neutron momentum reconstructed with LAC 
time-of- flight and the missing momentum for the reaction yd — > pX. 

To study the LAC response, the neutrons candidates were selected through 
the deuteron photo-disintegration 'yd — > pX imposing a cut in the missing 
mass distribution below the pion production threshold and the direction of 
the missing momentum oriented within the solid angle defined by the two 
LAC modules. The data were taken using the 4.0 GeV tagged photon beam 
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Fig. 11. The relative width of the difference between the neutron momentum re- 
constructed with LAC time-of-flight and the missing momentum for the reaction 
'yd — > pX as a function of P n . Empty dots: experimental data; triangles: Monte 
Carlo simulations. 



incident on a liquid deuterium target. At this energy the neutrons detected in 
the LAC modules come mainly from the exclusive reaction 'yd — *> pn. Typical 
neutron momenta vary between 0.7 GeV/c and 1.5 GeV/c. The comparison 
between the momentum P n reconstructed by the LAC time-of-flight and the 
missing momentum Px for neutron candidates is reported in Fig. 10. The 
distribution of the relative difference (P n — Px)/P n is correctly centered on 
zero. The width of this distribution has been studied as a function of P„ 
and compared to the Monte Carlo expectations. The results are reported in 
Fig. 11 where the empty dots represent the experimental data and the full dots 
the simulation results: the agreement is fairly good in the whole investigated 
momentum range. 



To evaluate the LAC efficiency we simply counted the fraction of candidates 
that generated a neutral hit in the LAC modules. In Fig. 12 we report the 
results (crosses) together with the Monte Carlo results (full dots) as a function 
of P n : as for the momentum resolution, there is complete agreement in the 
whole investigated momentum range. 
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Fig. 12. The neutron detection efficiency as a function of the neutron momentum. 
Empty dots: experimental data; Triangles: Monte Carlo simulations. 



6 Conclusions 



The two modules of the large angle calorimeter of the CLAS detector have 
been commissioned in 1998 and since then successfully operated. In this paper 
we have first discussed the calibration procedures proving that the detector 
geometry allows for simple energy and timing calibrating algorithms. These 
procedures are based on the detection of minimum ionizing particles (pions) 
that provide the relative energy and timing intercalibration. The absolute 
energy calibration for electromagnetic shower detection is obtained applying 
a correction factor as obtained from Monte Carlo simulations, and has been 
checked with the reconstructed ir° mass detecting one photon in the LAC 
module and one photon in the forward calorimeter. The dependence of the 
resolution from the deposited energy has been studied with protons completely 
stopped in the calorimeter. The timing absolute calibration is obtained setting 
(3 = 1 for the photons that originate from the ir° decay. The overall timing 
resolution for photons has been measured equal to 250 ps. Finally the response 
to neutron has been studied as a function of the neutron momentum. The 
obtained momentum resolution and the detection efficiency fairly agrees with 
the Monte Carlo simulations. In particular the neutron detection efficiency is 
larger that 20% for neutron momenta higher than 0.7 GeV/c. 

Therefore we can conclude that the two LAC modules show the expected 
performances. 
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